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It has been well known that quantum fluctuations induce a macroscopic phase transition from a 
superfluid to a Mott insulator phase driven by the repulsive potential energy in the ground state 
of dense bosonic systems on a lattice. We find a quantum phase transition from the homogeneous 
to an inhomogeneous Bose-condensate strongly affected or sometimes driven by the kinetic energy 
dispersion in dilute bosonic superfluids, which provides a clear identification of the superfluid state. 



The experimental realization of the quantum phase 
transition (QPT) from a superfluid liquid to a Mott in- 
sulator phase in an atomic gas trapped in an optical lat- 
tice [l| triggered remarkable experimental and theoreti- 
cal activity. This experiment heralded a new regime in 
exploration of the many body physics dominated by an 
interplay between atom-atom interactions and the boson 
kinetic energy as described by the Bose-Hubbard model 
. The competition between kinetic and interaction en- 
ergy terms in the underlying Hamiltonian has been con- 
sidered to be fundamental to quantum phase transitions 
not only in neutral but also in charged Bose liquids in the 
context of granular superconductors and Josephson junc- 
tion arrays iS] and preformed real-space electron pairs as 
bipolarons 

The Bose superfluid to Mott insulator transition takes 
place only at fixed density in a homogeneous dense sys- 
tem when the number of bosons is commensurate with 
the number of lattice sites. At first glance, one would not 
expect any phase transformation in a dilute Bose super- 
fluid well described by the Bogoljubov theory Sur- 
prisingly (as we shall show), this wisdom is not always ap- 
plicable to bosons in a periodic potential, whose ground 
state is described by a Gross-Pitaevskii-type (GP) equa- 
tion Ji, J,] including lattice, V{r), and interaction, U{r), 
potentials, 
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Here the condensate wave- functions V'(r) , which is also 
the order parameter, is normalized by the number of 
bosons as / (ir|-!/)(r)|^ = Nb, m is the boson mass, and /i 
is the chemical potential, which controls the total num- 
ber of particles in the system. One can solve equation 
([T|) variationally on a representative lattice with a large 
number of sites, N, using a complete set of orthogonal 
Wannier (site) functions w{r). Transforming the order 
parameter as ■!/'(r) = J2m '/'mw(r — m) we reduce ([J) to 
a discrete set of equations for the site amplitudes (/)m, 
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Here t(m) = J drw* {r)[-h'^V'^ /2m + V{r)]w{r - m) 



J J dvdv'U{v ~~ v')w*{v — m')u'*(r' - n')w{Y' — n)u'(r — m) 
is the matrix element of the interaction potential. 

We are interested in the dilute liquid regime far away 
from the Bose liquid-Mott transition. It is accustomed 
to keeping only the density-density interactions, C/™n" ~ 
^mn'5m,m"5n,n'- Solving the GP equation ^ is then 
equivalent to minimization of the energy functional. 
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Importantly, rescaling the order parameter as 0m — 
n^^^fm and the interaction as J7™n = Uji^/n yields a 
universal functional, i?(/m) = E{(f)„i)/n that is indepen- 
dent of n = Nt/N, the number of bosons per lattice 
valley: 



i(m - n) -I- ^(5„ 



1 



/n/n 



is the hopping integral, and 
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A solution to the tight-binding form of ([T]) corresponds 
to the minimum of this functional when the following 
constraint is imposed, l/mP = N. Solutions for dif- 
ferent particle densities are mapped on to each other by 
a simple rescaling of the interaction. 

The noninteracting part of Eq.(IT]) can be solved us- 
ing the Bloch eigenfunctions, and the Wannier states are 
constructed by summing the Bloch states with appropri- 
ate phase factors. Having obtained the Wannier states 
the hopping integrals and the interaction matrix elements 
can be evaluated for any periodic lattice potential V(r). 
When the lattice potential is deep enough, the nearest- 
neighbor hopping, ti dominates over the next-nearest 
(nn), t2 and the next-next (nnn), ^3 nearest neighbor tun- 
nehng Q. Lowering the lattice potential causes t2 and/or 
ts hopping integrals to be of the same order as ti. Also 
lowering the boson density reduces the rescaled repulsion 
between bosons. That allows us to investigate QPTs in 
the superfluid state which might be driven by the kinetic 
energy dispersion, i?(k) = i(m) exp(zk ■ m) rather 
than by the repulsion under certain conditions. 

We have randomly generated starting values for real 
site amplitudes, /m- The kinetic portion of the functional 
is calculated in the momentum space using the Fourier 
transformed values of the amplitudes at each lattice site. 
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FIG. 1: (a) Energy band dispersion and (b) inhomogeneous 
density when uoo ~ 4.0 in a deep square lattice potential 
with the hopping integrals ti — l,t2 — t3 = 0. c) Location 
of onset of inhomogeneity as t2 is varied. Gray indicates a 
homogeneous condensate, black inhomogeneous. Filled circles 
correspond to the locations of simulations run on 50 x 50 
lattices in the border region. 



Then it is inverse Fourier-transformed back into the Wan- 
nier space, where the effects of the interaction energy are 
apphed; note that the minimization procedure tends to 
fail in the absence of any repulsion. The functional (|4]) is 
minimized for a given value of fi using the NAG optimiza- 
tion routine E04DGF. The value of fi is then fine-tuned 
so that the constraint J2m l/mP = is imposed. There 
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FIG. 2: Energy band dispersion in a shallow square lattice 
potential with ti = l,t2 = 0, is = —0.4 (a), real space order 
parameter (b), and striped particle density (c) on a section of 
a 50 X 50 lattice. Here the weak on-site repulsion is woo = 0.02 
with no inter-site interactions. 



are virtually no significant size effects as verified by our 
simulations on = 25 x 25 and A^ = 50 x 50 lattices. 

Some of our results are illustrated in Figs. (|ll2l3p . 
Fig.(IT^) represents the single particle dispersion, E(k), 
for a deep square lattice potential, where ti is positive 
and t2 and ^3 are negligible. The level position in a single 
lattice valley is taken as zero, t(0) = 0, and all energies 
are measured in units of ii = 1. The order parameter and 
the density are uniform in this case, /m = 1, and there 
is no phase transition with increasing but moderate on- 



3 



site repulsion, Uqq < 1. Effects driven by the strength 
of the repulsion will be typically seen when uqq « 4.0 
or larger; these effects destroy the homogeneity of the 
condensate. Fig. ([TJd) shows an example of such an in- 
homogeneous condensate. Fig ^) shows how the value 
of the on-site repulsion uqo at which the homogeneity 
is destroyed decreases with the value of ^2- This corre- 
lates with the increasing shallowness of the kinetic en- 
ergy dispersion as t2 approaches —0.5, where the dis- 
persion becomes a pair of flat valley intercepting one 
another in a cruciform centered at the T point. This 
homogeneous-inhomogeneous transition occurs well be- 
low the critical value of uqo ~ 16 — 23 required for the 
superfluid-Mott insulator transition on a commensurate 
lattice [9I, [13, [m [12, [13] , and so is a property of the su- 
perfluid state. There are different irregular patterns in 
the inhomogeneous phase depending also on the intersite 
interaction uqi , which are reminicent but not identical to 
transformations of a homogeneous Bose condensate into 
a density wave superfluid [3, [3] with a tendency to the 
phase separation |16l |. 

Can one increase the magnitude of t2 further? In a 
conventional optical lattice, one cannot, as a result of 
the textbook theorem to the effect that the ground state 
wave-function of a single particle must not contain any 
nodes [13, [ill- This places hmits on the values of hopping 
parameters that are physically meaningful, since they 
cannot take values that give rise to additional minima in 
the kinetic dispersion of the lowest band of the Hamilto- 
nian. However, there are a number of feasible situations 
where the theorem is not applied due to many-body ef- 
fects [l^, [2^ or through the effects of internal degrees 
of freedom [2l|, [22)] . Recent experimental [1^, and 
theoretical work [25| has also suggested that condensa- 
tion may be induced in an excited but metastable p-band 
state. Also, this concern does not arise with composite 
bosons such as bipolarons [3|, since the single polaron 
band, where the condensation is occurring, is rarely the 
ground state band due to the Pauli exclusion principle. 
In this case the hopping parameters are determined by 
the symmetries of the wavefunctions of two components 
of the boson, so that there is no constraint on variation 
of any hopping parameter. 

With these considerations we turn to the results pre- 
sented in Figures [2] and [3] Fig. ([2]) represents a shal- 
low square lattice, where the magnitude of ts is com- 
parable with ti. At some value of the minimum of 
the single-particle band shifts from the F point with 
k = to 4 finite wave-vectors within the first Brillouin 
zone, ki = fc(l,l), ka = /c(-l,l), kg = fc(-l,-l), 
and k4 — fc(l, —1) where Bose condensation takes place. 
Fig. ((21) (the value of k is calculated below). The four min- 
ima are degenerate. A small repulsion removes the degen- 
eracy, so that the true condensate wave function in the 
Wannier space is one of the superpositions respecting the 
parity and time-reversal symmetry, oc cos(ki -m), 
/m oc cos(k2 • m), or /m oc cos(ki ■ m) ± cos(k2 • m), 
depending on the repulsion. In our example with the on- 
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FIG. 3: Energy band dispersion in a sliallow square lattice po- 
tential witli ti — ~l,t2 = —0.6 and fa = (a), the real space 
order parameter (b), and the checkerboard particle density 
(c) on a section of a 50 x 50 lattice. Here the weak on-site 
repulsion is uqo = 0.02 , and the nearest-neighbor interaction 
MOl = 0.01. 



site repulsion only the order parameter and the density 
of condensed bosons is striped along a diagonal direction, 
Fig. I©. Here the nnn tunneling plays the role of a pa- 
rameter driving QPT. One can find its critical value, t^c 
corresponding to the transition from the uniform to the 
striped condensate by calculating the effective mass m* 
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of a single boson on the lattice. With the positive ti = 1, 
zero ^2 and a negative ^3 the dispersion law is given by 
E{k) = -2[cos(fc2,) + cos(fcy)] - 2i3[cos(2fca;) + cos(2fcy)], 
where the lattice constant is taken as a = 1. Expanding 
in powers of k^^y one obtains i?(k) « — 4— 4i3 + /i^fc^/2r7i* 
with the effective mass m* — h^/{\ + 4t3), so that 
QPT appears at t^c = —0.25 when m* = cxd. The 
QPT order parameter is the wave-vector (or inverse pe- 
riod of the modulation), which above t^c is given by 
k = cos"-^ (1/41^3 1). It changes continuously from zero 
at the transition up to fc = 7r/2 at large \t^\, so that 
QPT is of the second order. The transition is driven en- 
tirely by the kinetic energy dispersion, rather than by the 
repulsion, so that it could be named a kinetic quantum 
phase transition (KQPT). A similar calculation for the 
case where ^2 is varied and held fixed at zero shows that 
there is a phase transition at t2c — —0.5 to a state with 
4 degenerate minima located at k = (tt, 0), (— tt, 0), (0, tt) 
and (0, — tt), and is discussed in relation to Fig. ([3]). Due 
to the discontinuous nature of this phase transition, it is 
first order. 

Quite generally, if the boson band dispersion has its 
minima at finite k the Bose condensate is nonuniform. 
Another example is the center-of-mass band dispersion 
of a small bipolaron composed of two holes on neigh- 
boring oxygen ions proposed as an explanation of the 
unusual symmetry and checkerboard modulations of the 
order parameter in cuprate superconductors [2l|. The 
nearest neighbor hopping integral is negative in this case 



due to the p-symmetry of oxygen orbitals, so that the 
minima of the bipolaronic band are found at the Bril- 
louin zone boundaries. Minimizing the energy functional 
([3]) with negative ti = —1 yields different patterns of the 
order parameter and the density depending on the inter- 
action matrix elements and longer range hopping terms. 
In this case we observe the relocation of 4 minima ini- 
tially located at k = (±7r, ±7r) to new positions within 
the first Brillouin zone, rather than a change from 1 to 4 
minima. An interesting example is the case of ^2 < —0.5 
and i3 — 0, where a weak nearest-neighbor repulsion uio 
can stabilize the d-wave checkerboard order parameter 
(see Fig. ([3])) - as anticipated in Ref.[2l| - provided that 
Moo is not too strong. Intriguingly, the same values of ^2, 
t3 and the repulsions also typically give a checkerboard 
state for the positive ti = 1. 

Cold atoms in optical lattices have provided an ex- 
cellent tool for investigating quantum phase transitions, 
but finding a reliable diagnostic criterion for superflu- 
idity is not straightforward [H, Our prediction of 
QPTs in dilute bosonic superfluids, which is strongly af- 
fected or even driven by the kinetic energy dispersion, 
opens up new perspectives on the unique diagnostic cri- 
teria of superfluidity. These states can most likley be 
detected through the use of optical probe spectroscopy 
and/or phase sensitive experiments. 
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